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ABSTRACT 


The application of aerodynamic theory for estimating the force and moments 
acting upon spinning airplanes is of interest. For example, strip theory has been 
used to generate estimates of the aerodynamic characteristics as a function of spin 
rate for wing-dominated configurations for angles of attack up to 90 degrees. This 
work, which had been limited to constant chord wings, is extended here to wings 
comprised of tapered segments. Comparison of the analytical predictions with rotary 
balance wind tunnel results shows that large discrepancies remain, particularly for 
those angles-of-attack greater than 40 degrees. 
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wing span 
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wing segment chord at the in-board edge 
lift coefficient 

rolling moment coefficient = qb area of segment 
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INTRODUCTION 


The use of parameter estimation in modeling aircraft dynamics has been quite 
successful for many mathematical models of flight. Parameter estimation is most 
readily applied when linear models representing small perturbations from straight 
equilibrium paths are appropriate. Flight data is most accurate in this regime and 
the mathematical model is the simplest. 1 

Parameter estimation becomes more complex in application to spinning aircraft. 
Modeling nonlinear aerodynamics , including rotational flow effects t is much more 

difficult and many more unknown parameters are introduced. 2 jn order to reduce the 

large number of unknowns it is helpful to apply strip theory of reference 3. Strip 
theory links the wing airfoil section characteristics to the rolling and yawing 
moment of the wing in spinning flight. "1 

In reference 1, strip theory provided a mathematical model that was used to 
determine the rolling moment of a wing in spinning flight. Calculated rolling moment 
forces due to the wing were about 50 percent larger than the experimental rotary 
balance spin-tunnel measurements of a wing-dominated aircraft. It is the purpose of 
this paper to expand the existing mathematical model of a spinning wing in order to 
more closely represent an aircraft in spinning flight, and to further explore the 
limitations and possibilities of the more general model. Specifically, the strip 
theory technique of reference 1 will be extended to wings comprised of tapered 
segments. The same limitation of reference 1 will be used in that the flow angle at 

each strip location is independent of the incremental lift at other locations. 


DISCUSSION 


In order to decrease the complexity of estimating the rolling moment due to 
spinning, the authors in reference 1 restricted their analysis to the rolling moment 
produced by an untapered wing of a wing-dominated aircraft. In this paper the 
approach is extended to wings of arbitrary planform by considereing a wing to be made 
up of sections of differing taper. 

Let us first consider the local flow characteristics for the general spanwise 
location y, shown in figure 1. 
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For wings having a constant taper, the wing chord can be represented by a linear 
equation: 

c = c - hy for y > 0 
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The equation for rolling moment for a single strip would be: 

dL = - | j(u - ry) 2 + (w + py) 2 j C N ^c o ± hy y dy^ 


For the entire wing the rolling moment becomes: 
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In order to easily represent aerodynamic data at high angles of attacks, the normal 
force coefficient (fig. 2) is given the form: 

C (a) = C m + C m sina 
N N o N sina 

It follows then that a single wing section over which the normal force equation 
is. applicable will have the following contribution to rolling moment: 
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After integrating, 
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The terras in the normal force equation, C», and C XI . , correspond to 

N o N sl na 

the local angle-of-attack ranges (see fig. 2) listed in Table 1 from reference 1: 
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The wing span locations having local angles of attack of -16, -10.5, 10.5 and 
16 degrees are determined by: 

w - u tan(a boundary) 

^boundary - p - r tan(a boundary) 

These will serve as limits of integration in the above equations if they fall in the 
confines of the panel being considered. If they do not, the boundaries of the panel 
will be used as limits. 

The program used to calculate the rolling moment of the wing using the above 
equations is listed in the appendix. It is a series of subroutines that will calcu- 
late the rolling moment coefficient of any tapered section of a flat wing given the 
following data: the two boundary chord lengths of each panel; the distance of these 
chords from the origin; the air density; the velocity of the aircraft; the wingspan; 
and the area of the wing. There is an option to calculate the rolling moment coeffi- 
cient of a single panel, or both symmetrical panels having the given dimensions. 

The spin subroutine accepts the dimensions of the panel and calculates the slope 
of the linear equation describing wing taper (h). It then computes the limits of 
integration along the panel. These limits are sent to the intermediate tests 
subroutine. Tests classifies the limits and sends only those that are within the 
bounds of the desired panel(s) to the panel subroutine. The panel subroutine does 
the actual rolling moment calculation of the panel between the limits using the above 
equations. The split subroutine is an optional subroutine which, given the 
dimensions of the wing, will split a wing into its component panels and send each 
panel in succession to the spin subroutine. 

With this program, a wing comprised of tapered panels can be modeled, panel by 
panel. Through a simple modification, the program can accumulate the total rolling 
moment of an aircraft wing due to each panel at a selected angle-of-attack. For the 
airplane shown in figure 3, this was done at an angle-of-attack of 14 degrees in 
order to obtain figure 4. Figure 4 is a plot of the total rolling moment coefficient 
of the wing of the aircraft, as well as the rolling moment coefficient of each of the 
wing's component panels as a function of nondimensional spin rate. The bottom curve 
of figure 4 represents the total rolling moment coefficient for the airplane of 
figure 3 at 14 degrees angle-of-attack. 
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The following flowchart is a diagram of the program: 



In a typical light, wing-dominated aircraft such as the one illustrated in 
figure 3, the panels that cause the greatest moment are the outer panels as is shown 
in figure 4 and in table 2. The upper curve in figure 4 represents the rolling 
moment contribution of the inner panels, the next curve represents the contribution 
of the middle panels and the third curve from the top is the contribution of the 

outer panels. This figure is for a fixed angle-of-attack while the rotation 
rate varies. On the other hand, table 1 shows the relationship between the panels 

when the rotation rate is fixed and the angle-of-attack is varied. The data of table 
1 and Figure 4 clearly show that the outer panels contribute from 78% to 97% of the 
total rolling moment. Of course, this is expected since these panels are larger than 
the others, have the longest moment arm, and experience the greatest variation in 

dynamic pressure. 

Figure 5 shows the improvement caused by taking into account wing taper as 
compared to the values obtained with a constant chord. There is significant 
improvement in the data, particularly at higher rates of rotation. The upper curves 
are the spin-tunnel test data. Obviously, improvements in the model must be made 
before the method can be considered acceptable. It is interesting to note (see 
figures 5 and 6) that there is little difference when the wing of the aircraft in 
figure 3 is simplified in the calculations to two large trapezoidal panels instead of 
six smaller ones. However, the multi-panel approach is more accurate and is 

applicable to the more general case. 
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In reference 1, it was noted that at angles-of-attack around 50 degrees the 
experimental rolling moments were autorotative at low rotation rates. The calculated 
data of reference 1 did not represent this phenomenon. The plot of 30 and 50 degrees 
angle-of-attack in figure 7 shows that the new calculated data does not show auto- 
rotative moments either. With the theory being used here, it would be impossible to 
obtain autorotative moments except over an angle-of-attack range of 10.5 to 16 
degrees since the slope of the line of normal force coefficient vs. angle-of-attack 
(fig. 2) is always positive except over this range. Note that figures 5 and 6 show 
an autorotative moment at low rates of rotation both in the test data and in the 
calculated data for 14 degrees angle-of-attack. However, an extension must be made 
to this simplified aerodynamic theory for higher angles-of-attack. 1 

The amount of error in the mathematical representation of a spinning wing has 
been decreased by describing the wing as a set of tapered panels. However, the 
errors are still large. The next step might be to consider the contribution of the 
tail section to the rolling moment. Since the program calculates the rolling moment 
of any tapered panel, the three tail panels could be input in order to determine the 
tail effects. The present method will compute the rolling moment for swept-wing 
configurations since rolling moment is independent of sweep. However, an extension 
of the model should also incorporate pitching moment. Of course, this method will 
not hold for aircraft where body effects cannot be neglected. The effects of the 
body would have to be considered by some other method such as the strip theory of 
reference 5. Improved estimates of aerodynamic moments would be expected if the 
induced flow effects on the flow angles were included in the formulation. Past 
results and future extentions promise further improvements in predicting the 
aerodynamic forces and moments of spinning airplanes. 


CONCLUDING REMARKS 


Mathematical representations of nonlinear phenomena such as the aerodynamics of 
a spinning aircraft are characterized by having large numbers of unknown parameters. 
Analytical methods such as strip theory can be used to reduce the number of unknown 
parameters. In this paper, strip theory is applied to compute aerodynamic forces for 
a wing composed of several variable taper trapezoidal panels in order to obtain a 
model structure which requires only the unknowns of the normal force equations. 
Although the error is decreased significantly by using strip theory in this manner as 
compared to approximating the wing as untapered, there is still much more to be done 
in order to analytically predict aerodynamic force of spinning aircraft. In order to 
extend the model further, many new parameters would have to be added. Also, it is 
clear that aerodynamic theory for angles-of-attack greater than 40 degrees must be 
improved since it is impossible to predict the results of spin-tunnel rotary balance 
tests with strip theory methods. 

Since the program that calculates the data is general enough to accept any wing 
panel of an aircraft, the revised model is currently useful in comparing the effects 
on a panel of changing parameters such as rotation rate, angle-of-attack, velocity, 
taper, etc. It is also useful for comparing aircraft components. However, the error 
between anlytical predictions and the experimental data is still too large to 
consider the strip theory representation to be an effective model of a spinning 
aircraft. 
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APPENDIX 


Listed in the following pages is the program used to calculate the rolling 
moment due to spin rate for an arbitrary planform wing. Inputs to the program are 
set in a short main program which calls the subroutines necessary for the 
calculations. 

The first page shows an example of the simplest case where a single panel is 
input to the program. Variable CHORD is the inboard chord length and CHTIP is the 
outboard chord length. D1 and D2 are spanwise distances from the center line of the 
fuselage to CHORD and CHTIP respectively. AREA refers to the area of the entire wing 
containing the panel, and SPAN is the wing span. RHO is air density and VEL is the 

velocity of the aircraft. The last integer tells the program whether to compute the 

rolling moment for panel with the given dimensions on the positive side of the 
aircraft (0), the negative side of the aircraft (1) or both (2). Normally, this 

value will be 2, except when isolation a single wing panel is desired. 

The second page shows a case where the panels of an entire wing will be input to 
the program. In this case, variable CHORD is the root chord, CHI and CH2 are the 
chord lengths at the point of wing taper change, and CHTIP is the chord length at 
the wing tip. D1 and D2 are the spanwise distances from the center line to CHI and 
CH2 respectively. FUSE is the width of the fuselage. 

As was mentioned in the text, the program can be modified to accumulate the 

total rolling moment of a multi-paneled wing. To do this, a one-dimensional array 
can be defined and placed inside the main loop of the spin subroutine such that each 
time spin is called with a new panel's dimensions, the nine values of the panel's 
rolling moment coefficient array (CLW) are added to the new array for a selected 
angle-of-attack (corresponding to an iteration of the main loop). For more than one 
angle-of-attack, a two-dimensional array would be necessary. 
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Figure I. Schematic Sketch ol a Spinning Wing 
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Figure 7. Iff act of rotation rata on rolling »o*ant coefficient for 30 and 
30 degree* angla of attack. 
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